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We make use of recent results in effective theory and higher-order perturbative calculations to im- 
prove the theoretical predictions of the QCD contribution to the top-quark pair production forward- 
backward asymmetry at the Tevatron. In particular, we supplement the fixed-order NLO calculation 
with higher-order corrections from soft gluon resummation at NNLL accuracy performed in two dif- 
ferent kinematic schemes, which allows us to make improved predictions for the asymmetry in the 
pp and tt rest frames as a function of the rapidity and invariant mass of the tt pair. Furthermore, 
we provide binned results which can be compared with the recent measurements of the forward- 
backward asymmetry in events with a large pair invariant mass or rapidity difference. Finally, we 
calculate at NLO+NNLL order the top-quark charge asymmetry at the LHC as a function of a lower 
rapidity cut-off for the top and antitop quarks. 

I. INTRODUCTION 

The forward-backward (FB) asymmetry in top-quark pair production in proton- antiproton collisions is an observable 
which originates from the difference in the production rates for top quarks in the forward and backward hemispheres 
[1, 2]. The total FB asymmetry was measured by the CDF and DO collaborations at the Tevatron [3-5]. The 
measurement can be carried out in the laboratory frame (pp frame) as well as in the center-of-mass frame of the 
top-quark pair (tt frame). The asymmetries in the two frames are defined as 

« _ N(y\ > 0) - N(yl < 0) 
FB N(y\ > 0) + N(yl < 0) ' 

where N is the number of events, i = pp (tt) indicates the laboratory frame (tt frame), and y\ is the top-quark rapidity 
in frame i. The measurements obtained by the CDF collaboration using 5.3 fb _1 of data are [5] 

Af B = (15.0 ± 5.5)% (pp frame) , 

A% = (15.8 ± 7.5)% (tt frame) . (2) 

The quoted uncertainties are derived from a combination of statistical and systematic errors. 

The production of top-quark pairs at hadron colliders is dominated by QCD. At the Tevatron, the charge conjugation 
invariance of the strong interaction implies that the difference in the production of top quarks in the forward and 
backward hemispheres is equivalent to the difference in the production of top and antitop quarks in the forward 
hemisphere. Therefore, in QCD the FB asymmetry is equivalent to the charge asymmetry. QCD predicts a non- 
vanishing contribution to the FB asymmetry starting at order a 3 in the squared amplitude. A contribution to the 
asymmetry arises if, in the interference of one-loop and tree-level diagrams, the top-quark fermionic line and the light- 
quark fcrmionic line are connected by three gluons [2] . The same is true in the case of the interference of two tree-level 
diagrams with three particles in the final state. The dominant effect is from the quark- antiquark annihilation channel, 
while a further, numerically small contribution to the asymmetry at order a 3 originates from the flavor excitation 
channel gq — > ttX, where X indicates additional partons in the final state. The gluon fusion channel does not 
contribute to the FB asymmetry at any order in perturbation theory, due to the fact that the gluon distribution is 
the same for protons and antiprotons. 

The total FB asymmetry predicted by QCD at the first non- vanishing order (which, for reasons discussed later in 
this paper, we will indicate as next-to-leading order (NLO)) is lower than the one measured at the Tevatron. Two 
recent evaluations using the formulas in [2] report the following values [6, 7], obtained using MSTW2008 NLO parton 
distribution functions (PDFs) [8], 

A% NLO = (4.8t°;|)% (PP frame) , 

^NLO = (7 4 +0.7 )% (ti - frame) , (3) 

The central values quoted in (3) refer to the choice [if = m t , where [if is the factorization scale and rat the top-quark 
mass, which is set to m t = 173.1 GeV. The errors originate from the scale variation in the range mt/2 < [if < 2mj. 
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These choices will also be adopted in the rest of the paper. Electroweak corrections enhance the prediction for 
the asymmetry by less than 10% of the central value [2, 9], and to separate uncertainties coming from electroweak 
calculations we do not include these corrections in our numerical results. As one can see, the discrepancy between 
the theory prediction and the experimental measurement in the pp frame is less than two standard deviations (2a), 
while in the tt frame the two values agree within ~ la, although the central value of the experimental measurement 
is a bit higher. 

In [5], the CDF collaboration measured the FB asymmetry in the tt frame as a function of the top-pair invariant 
mass M t f. After grouping the events in two bins corresponding to M t i < 450 GeV and M t i > 450 GeV, they found 
the asymmetry in the latter bin to be 

A$ B (M« > 450 GeV) = (47.5 ± 11.4)% , (4) 

which is more than 3a higher than the stated theoretical NLO prediction in [5] of (8.8 ± 1.3)% obtained using the 
MCFM program [10]. A measurement of the FB asymmetry in two bins of the rapidity difference yt — yi was also 
performed and again in that case the higher bin shows a tension with the NLO QCD prediction, although with larger 
experimental errors. Many attempts to explain these results in terms of new physics scenarios have been made, see e.g. 
[11-21]. The task is complicated by the fact that the new physics contributions should not spoil the good agreement 
between theory and measurements for the total pair-production cross section and the differential distribution in the 
pair invariant mass. 

Recently, we have calculated the top-pair invariant-mass distribution and the top-quark rapidity and transverse- 
momentum distributions in renormalization-group (RG) improved perturbation theory, which incorporates the re- 
summation of logarithmic corrections that become numerically large in the limit of soft-gluon emission [6, 7]. These 
calculations were carried out at next-to-next-to-leading-logarithmic (NNLL) accuracy. The resummed distributions 
were then matched to the NLO results to obtain predictions which have NLO+NNLL accuracy. Alternatively, the 
resummed results at NNLL order can be employed to generate approximate next-to-next-to leading (NNLO) pre- 
dictions [6, 7]. We note that working at approximate NNLO accuracy does not alter in a fundamental way any of 
the results obtained within the NLO+NNLL framework. On the other hand, we believe that in certain kinematic 
regions, for instance at large pair invariant mass, resummation is important. For these reasons, and for the sake 
of simplicity, in the following we focus on NLO+NNLL calculations. By integrating the differential distributions it 
is then straightforward to calculate the top-quark FB asymmetry both in the laboratory frame and in the tt frame. 
Starting from the double differential cross section in the pair invariant mass and scattering angle, we can also compute 
the M t i and rapidity-dependent asymmetries, which can be compared to experimental measurements. The main goal 
of this letter is to present the result of these calculations of the total and differential FB asymmetries in a systematic 
way. 

At the Large Hadron Collider (LHC), the symmetry of the pp initial state dictates that the rapidity distributions of 
the top and antitop quarks are symmetric and that the FB asymmetry vanishes. However, it was observed in [2] that 
at the LHC top quarks are preferably produced at larger rapidities than antitop quarks in the laboratory frame. Like 
the FB asymmetry at the Tevatron, this rapidity-dependent charge asymmetry is generated at order a 3 s in the squared 
amplitude, mainly through the asymmetric part of the quark- ant iquark annihilation channel. Therefore, potential 
new physics contributions would effect these two quantities in a correlated way, and the higher collider energy at the 
LHC gives it better access to distortions at higher rapidities. As a final application of our formalism, we evaluate 
at NLO+NNLL order the partially integrated charge asymmetry at the LHC, giving results as a function of a lower 
cut-off on the top and antitop rapidities in the laboratory frame. 



II. FB ASYMMETRY IN THE LABORATORY FRAME 



A. Total asymmetry 

The FB asymmetry in the laboratory frame can be calculated starting from the top-pair production cross section 
differential with respect to the top-quark transverse momentum px and rapidity y t . To do so, it is convenient to first 
define a total and differential asymmetric cross section via 
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where s is the square of the hadronic center-of-mass energy. To obtain the FB asymmetry in the laboratory frame 
one needs to calculate the ratio of the asymmetric cross section in (5) to the total cross section: 



app _ 



(7) 



In our phenomenological analysis we will consider two levels of perturbative precision for the asymmetric cross 
section and FB asymmetry. The first involves the differential cross section at NLO in fixed-order perturbation theory, 
the second the NLO calculation supplemented with soft-gluon resummation to NNLL order. In the laboratory frame, 
the resummed calculations are carried out using the IPIscet scheme introduced in [7]. Such predictions resum a 
class of logarithms which become large in the limit where S4 = (pt + k) 2 — m 2 is small (pf indicates the momentum 
of the antitop quark, and k the total momentum of additional partons in the final state), in which case real gluon 
emission is soft although the top pair is not produced at rest. While the phase-space integrals in (5) are in general 
sensitive to regions where S4 is not small, the threshold region is dynamically enhanced due to the sharp fall-off of the 
PDFs away from small values of S4. In fact, it was shown in [7] that the leading terms in the IPIscet threshold limit 
reproduce essentially the entire NLO correction to the differential cross section in the quark channel, which implies that 
rcsumming the corrections which become large in this limit is an improvement on the fixed-order calculation. While 
in the fixed-order counting the asymmetric cross section first arises from the NLO calculation of the differential cross 
section, in the resummed counting it first appears at NLL order. The NLO+NNLL calculation is thus a refinement 
on the leading term, and will be considered our best prediction. 

Before illustrating our results for the FB asymmetry, we need to clarify an important point concerning our convention 
for counting orders in the perturbative expansion. We calculate Afb itself as a perturbative expansion in a s , using a 
fixed-order or logarithmic counting as appropriate. For example, the first non- vanishing contribution to the asymmetry 
in fixed-order perturbation theory is obtained by calculating the numerator in (7) to order a\ and the denominator 
to order a 2 . The resulting asymmetry is of order a s , which we will refer to as NLO, with reference to the order 
at which the differential distributions in (5) are calculated relative to a 2 . Similarly, in RG-improved perturbation 
theory, the first non- vanishing contribution to A?b is obtained by calculating both the numerator and the denominator 
at next-to-leading-logarithmic (NLL) order; the resulting asymmetry is then counted as NLL. This convention was 
already adopted in [6, 7]. There arc two counting schemes in the literature (in fixed-order perturbation theory) which 
arc different from ours. The first one also treats Aps itself as a perturbative expansion, but counts the order a s 
contribution as LO, and so on, as in [I, 2]. The second one treats the numerator and the denominator as separate 
perturbative series and does not further expand the ratio, as adopted in the quoted MCFM results in [5]. We note 
that while the first scheme differs from ours only by name, the second scheme leads to different numerical results. As 
discussed in [6], the NLO+NNLL results are considerably more stable with respect to the choice of scheme than the 
NLO results. 
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TABLE I: The asymmetric cross section and FB asymmetry in the pp frame. The first error refers to perturbative uncertainties 
estimated through scale variations as explained in the text, and the second error in the MSTW2008 case is the PDF uncertainty. 



Our results for the total FB asymmetry in the lab frame are shown in Table I. As explained above, to calculate 
each entry in the table the numerator and the denominator in (7) are evaluated at the order indicated in the leftmost 
column, and then the ratio itself is expanded in powers of a s up to the appropriate order. The central values are 
obtained by fixing the factorization scale at /x/ = m t , and the scale uncertainties are estimated by varying p f between 
m t /2 and 2m t - 1 In the resummed calculations also the hard and soft scales are varied according to the procedure 



1 Although wc use mt = 173.1 GeV throughout the analysis, the asymmetry is rather stable under the exact choice of mt. For instance, 
at mt = 160 GeV the default value for the asymmetric cross section at NLO with MSTW2008 PDFs changes to Act FFs = 0.384 pb, but 
the asymmetry itself changes only by a small amount to Apg = 4.67%. 
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described in [7], and the effect of these variations is included in the scale uncertainty shown in the tables. In the 
first column of the table we use MSTW2008 PDFs and estimate the PDF uncertainties by iterating through the 90% 
confidence level (CL) sets. Here and throughout this work we use NLO PDFs for the NLO prediction, and switch to 
NNLO PDFs for the NLO+NNLL calculation. The PDF uncertainties for the asymmetry, expressed as a percentage 
of the central values, are about half as large as those for the asymmetric cross section. This is due to cancellations 
in the ratio. In addition to results with MSTW2008 PDFs, we also show those obtained using CTEQ6.6 [22] and 
NNPDF2.1 [23] PDFs. In those cases the PDFs are based on a NLO fit so that the same set is used in both the NLO 
and the NLO+NNLL calculations. We note that the results for the asymmetry obtained with the different PDF sets 
are well within the PDF uncertainties estimated through the MSTW2008 results. 

Adding soft-gluon resummation at NNLL accuracy produces results for the asymmetric cross section which are 
numerically consistent with the NLO results for fif — m t , while the scale uncertainty is reduced by more than a factor 
of 2. The central value for the FB asymmetry does not change significantly with respect to the NLO predictions, and 
also in this case the scale uncertainties are reduced. We can therefore conclude that the discrepancy between theory 
and experiment cannot be explained with the effect of higher-order QCD corrections on the theory side, at least not 
those related to soft-gluon resummation. 

The NLO+NNLL calculation represents the most accurate determination of the QCD contribution to the asymmetry 
that can be obtained at present. However, it is important to keep in mind the uncertainties related to yet higher-order 
corrections and how they could be reduced. Concerning soft gluon resummation, the most important improvement 
beyond NNLL accuracy would be the calculation of the hard and soft matching functions at NNLO order. This requires 
the soft plus virtual corrections at NNLO in QCD, and would fix the (5-function contribution to the differential cross 
section at that order. In that case, the further improvement of matching the resummed results to fixed order would 
require power-suppressed effects related to hard gluon emission. The arguments based on the dynamical enhancement 
of the threshold region and the confirmation of this mechanism through the numerical results at NLO, imply that 
power corrections to the soft limit are small, so we expect the more important effect to be the calculation of the soft 
plus virtual corrections. We estimate uncertainties related to both types of corrections through the standard method 
of scale variations and the numerical results indicate that these higher-order effects are moderate. However, this 
statement can of course never be certain without the actual calculation of the higher-order pieces. Similar comments 
apply to all other quantities obtained in this work. 



B. Rapidity-dependent asymmetry 



As experimental measurements become more precise, differential quantities such as the rapidity-dependent asym- 
metry can be compared with theoretical predictions. Using quantities defined in (5), we can write this differential 
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In Figure 1 we show results for the differential asymmetric cross section and the FB asymmetry as functions of the 
rapidity at NLO and NLO+NNLL order, using MSTW2008 PDFs. Here and below, the differential cross sections at 
NLO are obtained using a private NLO version of MadGraph [24] . The bands refer to uncertainties associated with 
scale variations. One observes that the NLO+NNLL band for the asymmetric differential cross section displayed in the 
left-hand panel is contained within the NLO band over the entire range of y t values shown in the figure, and that the 
scale uncertainty of the NLO+NNLL asymmetric cross section is smaller than the scale uncertainty obtained in the 
NLO calculation. On the other hand, the differences between the error bands in the NLO and NLO+NNLL results for 
the FB asymmetry shown in the right-hand panel are much smaller due to cancellations in the ratio which make the 
NLO FB asymmetry considerably more stable than the corresponding asymmetric cross section. We will encounter 
this feature repeatedly in the tt-frame calculations which follow. One observes that the form of the rapidity-dependent 
asymmetry, which is an increasing function with respect to yt, is very stable under higher-order corrections. 



III. FB ASYMMETRY IN THE tt FRAME 



A. Total asymmetry 



Our studies of the FB asymmetry in the tt rest frame will use as a fundamental quantity the top-pair production 
cross section differential with respect to the pair invariant mass and the top-quark scattering angle in that frame. We 
thus define an asymmetric cross section as 
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and the total FB asymmetry in the tt frame is then given by 

FB — ~ 



(9) 



(10) 



As in the previous section, we will study the FB asymmetry at both NLO and NLO+NNLL accuracy. In the tt 
frame, the resummed calculations arc carried out using the PIMscet scheme introduced in [6]. The differential cross 
section calculated in that work includes the resummation of partonic threshold logarithms up to NNLL order directly 
in momentum space, and thus extends the previous results from [25], which included resummation to NLO+NLL 
accuracy in moment space. Such predictions resum a tower of logarithms which become large when the square of the 
pair invariant mass is close to the partonic ccntcr-of-mass energy: z ~ M^/s — > 1. In this limit real gluon emission 
is soft although, as in the case of 1PI kinematics, the top quarks are not necessarily produced at rest. While the soft 
limit z — > 1 can be enforced kinematically by the restriction to large pair invariant mass, the phase-space integrals 
in (9) are in general sensitive to regions where this is not the case. However, the analysis of [6] showed that events 
near the partonic threshold provide the numerically dominant contributions to the differential cross section, even at 
low values of the invariant mass near the peak of the distribution, due to dynamical threshold enhancement from the 
PDFs [26], which indicates that resumming threshold logarithms is an improvement on the fixed-order calculation. 
The suppression of power corrections to the soft limit is also backed up by the negligible numerical difference between 
the calculation of the asymmetry in the tt frame and the partonic center-of-mass frame. The results in these two 
frames coincide in the threshold limit z — > 1, and differ by only about 1% at NLO [6], due to very small corrections 
from hard gluon emission. 2 



2 In [6], results at NLO+NNLL order for the total asymmetry were actually given in the partonic center-of-mass frame. In the present 
work, we have eliminated this small mismatch with the experimental measurements by working in the it frame. 
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TABLE II: The asymmetric cross section and FB asymmetry in the tt rest frame. The first error refers to perturbative 
uncertainties estimated through scale variations, and the second error in the MSTW2008 case is the PDF uncertainty. 



Our numerical results for the total asymmetric cross section and FB asymmetry are summarized in Table II. As 
was the case in the laboratory frame, the scale uncertainties in the asymmetric cross section are roughly halved at 
NLO+NNLL order compared to NLO. The scale uncertainties in the FB asymmetry, on the other hand, actually 
increase slightly after adding the resummation, while the central values are nearly unchanged. We note however that 
the resummcd results are more stable with respect to the scheme for expanding the ratio defining the FB asymmetry [6] , 
and in that sense are more reliable than the NLO predictions. Moreover, the resummed results for both the asymmetric 
and total cross sections are more stable under scale variations than their fixed-order counterparts. One should therefore 
be cautious of the rather small scale uncertainties in the NLO calculation of the FB asymmetry, which result from 
large cancellations in the ratio not observed in the resummed result. We again show the PDF uncertainties using the 
MSTW2008 PDFs at 90% CL, and the central values and scale uncertainties from the CTEQ6.6 and NNPDF2.1 sets. 
All comments from the previous section concerning the reduction of PDF errors in the FB asymmetry compared to 
the asymmetric cross section, and the good agreement between the different PDF sets, are also true in this case. 



B. Invariant-mass dependent asymmetry 

As mentioned in the introduction, the recent measurement of the asymmetry at high values of the pair invariant 
mass shows a large deviation from the value predicted by QCD at NLO. In order to examine the effects of soft-gluon 
resummation on this observable, we first extract from (9) the invariant-mass dependent asymmetry as 

( da \ ( da \ dAa£ B 
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Results for A§ B (Mu) and dAap B /dMu at NLO and NLO+NNLL order are shown in Figure 2, where the bands reflect 
uncertainties originating from scale variations. The figure shows that the asymmetry increases with the invariant mass 
and can reach nearly 40% at M tl = 1200 GeV. These results are obtained with the default MSTW2008 PDFs and do 
not include PDF uncertainties. An analysis shows that the relative PDF error for dAa^ B /dMu increases slightly with 
increasing M tf , from 7% (400 GeV) to 9% (1200 GeV). In contrast, the relative PDF error for A$ B {M tl ) is rather 
small and even decreases with M ti , from around 2% (400 GeV) to 1% (1200 GeV). 

It is a well-known fact that electroweak corrections start to become more important for the differential cross section 
at high pair invariant mass, due to the presence of Sudakov logarithms. For instance, at M t i ~ 1 TcV, the electroweak 
corrections to the differential distribution at the Tevatron are roughly at the -5% level [27, 28], even though for the 
total cross section they are negligible. On the other hand, the electroweak corrections to the asymmetric cross section 
given in [2] do not contain Sudakov logarithms, and an estimate shows that their distortion of the QCD contribution 
is roughly independent of the invariant mass. We therefore do not expect the electroweak corrections to significantly 
alter our results for the FB asymmetry even at high values of pair invariant mass. 

We hope that better statistics can eventually lead to a detailed comparison of experimental results with the asym- 
metry curve in Figure 2. At present this is not possible, but the CDF collaboration has measured the invariant-mass 
dependent asymmetry by separating the events into a high invariant-mass bin (M t j > 450 GeV) and a low invariant- 
mass bin [Mu < 450 GeV) [5]. From Figure 2, one can see that this choice roughly divides the total asymmetric 
cross section equally between the two bins, although most of the asymmetric cross section in the high invariant-mass 
bin originates from the region close to 450 GeV. More precisely, the region 450 — 600 (700) GeV captures more than 
75 (90)% of the total asymmetric cross section in the high invariant-mass bin. To compare with the CDF results, we 
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FIG. 2: Left: The asymmetric cross section dAa^ B /dM t t as a function of the invariant mass at NLO and NLO+NNLL order. 
Right: The asymmetry Ap B (M t j). The bands show the uncertainties related to scale variation as explained in the text. 
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FIG. 3: The asymmetry in the high and low invariant-mass region as measured in [5], compared to our predictions at 
NLO+NNLL order. The bands in the NLO+NNLL results are related to uncertainties from scale variation, while the NLO 
result in the higher bin is evaluated at fif — m t . 



evaluate the binned asymmetry 
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for M t j < 450 GeV and for M t j > 450 GeV. Our findings are given in Table III, along with their visual representation 
in Figure 3, which shows the NLO+NNLL calculation with an error band from scale variations along with the default 
NLO number in the high invariant-mass bin. In both bins, the NLO+NNLL predictions for the asymmetric cross 
sections have considerably smaller scale uncertainties than the NLO ones, but the results for the FB asymmetries are 
essentially unchanged. As with all other results obtained in the ti frame, the scale uncertainties in the FB asymmetries 
are larger in the NLO+NNLL calculation that at NLO. However, if we had not expanded the ratio, the predicted FB 
asymmetry in the high invariant-mass bin would be 9.0% at NLO and 10.6% at NLO+NNLL order 3 , showing the 



3 Using MSTW2008 PDFs as an example. 
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FIG. 4: Left: The asymmetric differential cross section dAap B /dAy. Right: The asymmetry Ap B (Ay). The bands show the 
errors related to scale variation as explained in the text. 



stability of the resummed results under this change of systematics. 
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TABLE III: Comparison of the low- and high- mass asymmetry Ap B with CDF data [5], along with results for the asymmetric 
cross section. The errors in the QCD predictions refer to perturbative uncertainties related to scale variation. 



We now turn to a discussion of the PDF uncertainties in the binned results, in this case deviating slightly from our 
usual procedure. The reason is that to compute the PDF uncertainties for the binned asymmetry at NLO+NNLL 
order, we need to run the Monte Carlo program MadGraph at NLO for each of the PDFs in the error set, which is 
rather time consuming. As a compromise, we have estimated the PDF uncertainties using only the pieces of the NLO 
calculation which are leading in the threshold limit M?t/s — > 1, as obtained in [6]. Since these leading pieces alone 
account for the bulk of the NLO FB asymmetry, the relative PDF uncertainties obtained from these terms should 
provide a good approximation to those in the full NLO and NLO+NNLL results. For the MSTW2008 set, we find a 
relative PDF uncertainty of about 7% for the asymmetric cross section and about 2% for the FB asymmetry at 90% 
CL, in both the low and high invariant-mass bins. 

Our calculations show that neither higher-order corrections from soft-gluon resummation nor the inclusion of a 
systematic uncertainty coming from PDF usage reduces in any significant way the current discrepancy between 
theory and experiment for the FB asymmetry in the high invariant-mass bin, which remains above the 3tr level when 
using our NLO+NNLL calculations. 



C. Rapidity-dependent asymmetry 

A further observable of interest is the rapidity dependence of the FB asymmetry in the ti frame. In practice, 
experiments measure the asymmetry as a function of the pair rapidity difference Ay — y t — yi [5]. We can calculate 
the differential cross section in this variable from the results in PIM kinematics by using that, up to power corrections 
which vanish in the soft limit, 



1 + cos 



1y/l - Ami/Ml 



Ay = In ( V , I ■ (13) 
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1 

Ay 

FIG. 5: The asymmetry for Ay < 1 and Ay > 1 as measured in [5], compared to our predictions at NLO+NNLL order. The 
bands in the NLO+NNLL results are related to uncertainties from scale variation, while the NLO result in the higher bin is 
evaluated at fj,f — m t . 
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TABLE IV: Comparison of Ap B for Ay < 1 and Ay > 1 with CDF data [5], along with the asymmetric cross section. The 
errors in the QCD predictions refer to the uncertainties related to scale variation. 

After changing variables from the scattering angle to the pair rapidity difference, we express the asymmetric cross 
section as 



a4b 



dAy 



dAy 



./,,, dMtl dM ti dAy 

( da \ ( da \ 
\dAy~) F ~ \dAy~) B 



d 2 a pp^ttX 

M mi„ dM ti dAy 



Ay=-Ay 



dAy 



dAa tl 



FB 



(14) 



where 



Ay 4 



In 



1 + y/l - 4m t 2 /s 



and M£ in = 2m t cosh (Ay/2) . 



1 - y/l - 4m 2 t / S/ 

Using these definitions, we can also introduce the Ay-dependent asymmetry 



(15) 



'FB 



Att 

il 



FB 



(Ay) = 



dAy 



da \ f da \ 
'dAy') F + \dAy) B 



(16) 



and binned asymmetries analogous to (12), where the numerator and denominator of the above expression are inte- 
grated over a range in Ay. Note that the integration region above implies that higher values of Ay correspond to 
higher values of M t j. For example, the restriction Ay > 1 used in the binned analysis below corresponds to events 
with M t i > 390 GeV. 

We show results related to the rapidity dependence of the FB asymmetry in Figures 4 and 5, and in Table IV. In 
all cases we use MSTW2008 PDFs. The more detailed results in the Figure 4 show the differential asymmetric cross 
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section along with the pair-rapidity dependent FB asymmetry, in the form of bands related to uncertainties from scale 
variations. While resummation stabilizes the asymmetric cross section compared to NLO, there is little effect on the 
FB asymmetry. The results for the binned asymmetry are given in Table IV, along with their visual representation in 
Figure 5, which shows the NLO+NNLL result with an error band from scale variations along with the default NLO 
number in the higher bin. For events where Ay < 1, the QCD prediction is in agreement with the CDF measurement 
[5]. In the bin where Ay > 1, the predicted asymmetry is lower than the measured one by ~ 1.5er. Again in this case, 
soft-gluon resummation changes the NLO predictions only slightly. 



IV. CHARGE ASYMMETRY AT THE LHC 



The Tevatron results for the FB asymmetry at high pair invariant mass and rapidity hint at a discrepancy with the 
Standard Model. It would of course be desirable to study the physics responsible for this effect through measurements 
at the LHC. 

The total and differential FB asymmetries at the LHC vanish, because of the symmetric initial state. However, 
while charge conjugation invariance of the strong interaction implies that the rate for the forward production of top- 
quarks is equal to the rate for backward production of antitop quarks at the Tevatron, this is not the case at the LHC. 
At a proton-proton collider the total rate for top and antitop production in the forward or backward hemisphere is 
equal, but at a given rapidity the rates differ. In fact, at large (small) rapidities the rate for top-quark production is 
noticeably larger (smaller) than that for antitop production [2], so although there is no FB asymmetry at the LHC 
there is a differential charge asymmetry. Like the FB asymmetry at the Tevatron, this charge asymmetry at the LHC 
is related to the asymmetric part of the qq partonic cross section, implying a direct correlation between potential 
new physics contributions to the two measurements. The charge asymmetry at the LHC is generally smaller than 
the FB asymmetry at the Tevatron due to large contributions from the gg channel to the charge-symmetric part of 
the differential cross section, but the rapidity reach at the LHC is larger and the charge asymmetry thus provides 
complementary information. 

In this section we study the simplest realization of a charge asymmetry at the LHC, namely the rapidity-dependent 
quantity in the laboratory (pp) frame. In particular, we focus on the partially integrated charge asymmetry and 
charge-asymmetric cross section, where we impose the restriction y > y cut on the differential cross section. We define 
these through 



A p c p (y C ut) = 



,7.„ t * \ d Vt dyt 



Vt=Vtj 



cut } 



y i I da pp ^ tx t d<r pp ^ lXt 

d Vt . , 

dy t 



.'/cut 



dyi 



F* dy ( dcjPP ^ tXt 
Jycut Vt \ dy t 



d G PV^iX t 



dyt 



(17) 
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with yf as in (6). 

We can study the partially integrated charge asymmetry and asymmetric cross section (17) at NLO and NLO+NNLL 
order using the results for the differential cross section in 1PI kinematics obtained in [7]. The results generated with 
MSTW2008 PDFs are shown in Figure 6, where the bands reflect uncertainties related to scale variations carried out 
with the same procedure as at the Tevatron. For the charge-asymmetric cross section shown in the left-hand panel of 
the figure, the main effect of the resummation is to decrease the scale dependence of the result to a relatively small 
region of the NLO error band. For the partially integrated asymmetry shown in the right-hand panel of the figure, 
the resummation is a mild effect up to y cut ~ 1.5, but substantially reduces the asymmetry at higher values of the 
cut. However, due to large K factors in the gluon channel at the LHC, the uncertainty band for the NLO curve is 
very sensitive to whether one consistently expands the asymmetry ratio in (17). If we had deviated from our normal 
procedure and had not expanded the ratio, instead evaluating the denominator at NLO order, the NLO band would 
actually overlap quite well with the NLO+NNLL results shown in figure. The NLO+NNLL result is largely insensitive 
to this change of systematics-the result where both the numerator and denominator are evaluated at NLO+NNLL 
order is within the error band shown in the figure. 

The partially integrated charge asymmetry vanishes for y cut = 0, and becomes progressively larger at higher values 
of the cut. However, the charge asymmetric cross section shown in the left-hand side of the figure is very small 
at higher rapidity values and the experimental measurement is difficult. A reasonable way to compare theory and 
experiment in this case would be to perform a measurement in a high-rapidity bin with y > y cut ~ 1-1.5. In such a 
bin the Standard Model charge asymmetry is predicted to be only slightly different from zero, less than 2% depending 
on the exact choice of the cut, so any appreciable charge asymmetry in an experimental measurement would be a 
clear sign of a new physics contribution to the high-rapidity region of the distribution, which is already hinted at by 
the Tevatron measurements. 

Note that other partially integrated differential distributions can be used to study the phase-space dependence of 
the charge asymmetry. For instance, one can consider the so-called central charge asymmetry introduced in [29], 
which imposes a cut \y\ < yc along with a restriction to relatively high pair invariant mass in order to reduce the 
contribution from the gg channel to the symmetric part of the partially integrated cross section. Unfortunately, it 
is beyond the scope of the paper to perform an NLO+NNLL analysis for such an observable, because it requires a 
combination of PIM and 1PI kinematics which is not possible to derive starting directly from the results of [6, 7]. 

V. CONCLUSIONS 

The total top-quark FB asymmetry measured at the Tevatron is not in good agreement with the predictions 
obtained at the first non-vanishing order in perturbative QCD. The tension between theory and experiment is about 
two standard deviations for the asymmetry measurement in the laboratory frame, and approximately la in the top-pair 
rest frame. Moreover, the measurement of the asymmetry in the large pair invariant-mass region (M t j > 450 GeV) by 
the CDF collaboration [5] differs from the leading-order prediction in QCD by more that three standard deviations, 
and there is also a discrepancy between theory and experiment for the FB asymmetry at values of the top-pair rapidity 
difference Ay > 1. 

The calculations of the double differential distribution for the top-quark pair production in two different kinematic 
schemes, carried out in [6, 7], allow us to evaluate the FB asymmetry in both the laboratory frame and in the top-pair 
rest frame. The predictions obtained in this way include the full NLO corrections to the total and asymmetric cross 
sections, whose ratio determines the FB asymmetry, as well as the resummation of soft-gluon emission effects up to 
NNLL accuracy. Such NLO+NNLL results represent the most accurate determination of the QCD contribution to the 
asymmetry that can be obtained to date. Studies of the NLO results indicated that power corrections to our results 
obtained in the soft limit within the effective theory approach are small, due to the mechanism of dynamical threshold 
enhancement. The full NNLO calculation would be very useful in quantifying this at the next order in perturbation 
theory, but we do not expect such power corrections to change our predictions in any significant way. 

Tables I and II show that the numerical impact of the NNLL corrections on the total asymmetry is modest, less than 
3% of the central value in all cases considered. It must be observed that while the total and asymmetric cross sections 
calculated at NLO+NNLL order have smaller scale uncertainties than the corresponding quantities calculated at 
NLO, the same is not always true for the FB asymmetry, indicating that the errors at NLO are affected by accidental 
cancellations and do not reflect the true perturbative uncertainties. After the inclusion of the NNLL corrections, the 
predicted values of the asymmetry remain lower than the measured ones by <~ 2a in the pp frame and ~ la in the tt 
frame. We have also shown that PDF uncertainties are much reduced for the FB asymmetries compared to those for 
the asymmetric or total cross sections. For instance, the total asymmetry calculated with MSTW2008 PDFs differs 
from the one obtained from CTEQ6.6 or NNPDF2.1 PDFs by only a few percent. 

The NLO+NNLL calculation of the asymmetry in the high invariant-mass bin M t i > 450 GeV increases the NLO 
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prediction only by a very small amount, about 3% of the central value. The perturbative and PDF uncertainties of the 
calculation are much smaller than the experimental errors. Hence, if the Tevatron measurement of the FB asymmetry 
in the high invariant-mass bin is confirmed as the experimental error is reduced, it will become a compelling signal 
for physics beyond the Standard Model. 

Given the situation at the Tevatron, it is well motivated to measure observables sensitive to the top quark asymmetry 
also at the LHC. Since the FB asymmetry vanishes at the LHC, one of the possible choices consists in measuring the 
charge asymmetry partially integrated in a given rapidity region. Starting from the NLO+NNLL resummed rapidity 
distributions for the top and antitop quarks, we calculated the charge asymmetry as a function of a lower rapidity 
cut in the laboratory frame. A careful choice of the lower bound on the rapidity integration region can allow one 
to use this observable to test for the presence of new physics in the top-quark sector at the LHC, thus providing 
complementary information to the Tevatron measurements. 
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